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Abstract. We study the bipartite entanglement between a sub-system of size I and 
the rest of the system of total size L as it occurs in a spin-1 AKLT chain subject to open 
boundary conditions. In this case the ground-state manifold is four-fold degenerate and 
there is strong dependence on the parity of the number of spins, L. We present exact 
analytical results for the von Neumann entanglement entropy, as a function of both 
the size of the sub-system, I, and the total system size, L, for all four degenerate 
ground-states for both odd an even L. In the large /, L limit the entanglement entropy 
approaches ln(2) for the = ±1 while it approaches twice that value, 21n(2), for the 
Sj. = states. In all cases, it is found that this constant is approached exponentially 
fast defining a length scale £ = l/ln(3) equal to the known bulk correlation length. 
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1. Introduction 

Quantum entanglement as it occurs in quantum spin chains is a property that has 
recently been under intense study [U [21 El El E|. The entanglement of a system can 
provide information about the properties of that system, and long-distance entanglement 
is thought to be necessary for applications such as quantum teleportation [El [7] and 
quantum cryptography. [8] Exact results for the bipartite entanglement is from this 
perspective of considerable value. 

One system displaying entanglement is that of an S — 1 antiferromagnetic chain. 
[9] One generalized Hamiltonian for such a chain is given by: 

H = JJ2lS t ■ S i+1 - P(Si ■ S i+1 ) 2 }, (1) 
1=1 

where is a dimensionless parameter describing the biquadratic coupling. When 
= —1/3, the system is at the Affleck-Kennedy-Lieb-Tasaki (AKLT) point, where 
the ground state of the system corresponds to a system where each S = 1 spin 
is represented as two S = 1/2 spins, and S = 1/2 spins from neighboring sites 
are combined into a singlet. [101 [11] One and two-site entanglement at the AKLT 
point and for more generalized models has been studied extensively with periodic 
boundary conditions. p21 [131 E] When subject to periodic boundary conditions the 
hamiltonian, Eq. ([1]), has a nondegenerate singlet ground state at the AKLT point. [15] 
Some measures of entanglement has also been studied for the case of open boundary 
conditions [HI [TTl HH HH [201 [SU [221 E3] , in this case the ground state of the system is 
four- fold degenerate [TQl LTl] , consisting of a singlet state, S = 0, as well as a triplet state 
with S = 1, S T = 0, ±1. An interesting quantity to study is the bipartite entanglement 
entropy S(l,L), the von Neumann entanglement of a subsystem of the chain with the 
rest of the chain: 

S(l,L) = -Trp\ogp, (2) 

where p is the reduced density matrix for the subsystem of size I within the total system 
of length L. Similar calculations have also been performed for S = \ systems. [2^ l l25 | [26] 
In physical systems well characterized as S = 1 spin chains such as NENP [27J and 
Y 2 BaNi05 [2H] the biquadratic term is negligible, (3 = 0, and impurities likely cut the 
chains thereby effectively imposing open boundary conditions and restricting the length 
of such finite chain segments. The presence of the open boundaries has the peculiar effect 
of inducing 5=1/2 excitations localized at the ends of the chain segment [29]. The 
physically most relevant point, (3 = 0, is in the same phase as the AKLT point, the so 
called Haldane phase. Within the Haldane phase, for (3 ^ —1/3, the four- fold ground- 
state degeneracy is lifted for finite L and is replaced by an exponentially low-lying 
triplet state above the singlet ground-state when the length of the system is even. For 
odd length systems the picture is reversed and the triplet state is lowest. A complete 
characterization of the entanglement as it occurs for all 4 states in the ground-state 
manifold would therefore be of interest. While S(l, L) has been studied at the AKLT 
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point for periodic boundary conditions by Hirano et al. [H] the only result, by Alipour 
et al. [23] , for the physically more interesting case of open boundary conditions is for the 
special case / = 1 with Sj< = ±1. In the following, we present analytical results for the 
bipartite entanglement entropy for the AKLT system with open boundary conditions, 
as a function of both the size of the total system, L, and the size of the subsystem, I. 
We explicitly present results for all four states for both even and odd length systems. 



2. Calculations 



In order to facilitate the calculations it is convenient to write the ground-state wave 
function of Eq. (CE]) in the following manner: [THl ELT1 EQj EU [32] 

In the above equation, |0)j, and |— )j are the states of the S=l spin at site i. The 
matrix gives the four ground states of the system. The upper right and lower left entries 
correspond to the states with magnetization S = 1, — 1 and S — 1, — — 1. The 
two Sj. = are as written not part of total spin multiplets and in light of the splitting 
of the ground-state manifold away from the AKLT point it is therefore of interest to 
form total spin eigenstates which is conveniently done by defining: 

\y >triplet=^(lp-Sl(<lP)) [ ' 

where ip is a diagonal entry in the wave function matrix, and SI is the spin inversion 
operator. We can use this wave function to calculate S(l,L). First we break the spin 
chain up into two subchains, A and B. We denote the number of spins in A by I, 
compared to the total number of spins L. To find S(l, L), we must first find the reduced 
density matrix, given by 

I' -j X/' '/'<'./• ( 5 ) 

j 

where i and i' run over all possible configurations of the subsystem A, and j runs over 
all possible configurations of the subsystem B. We can then diagonalize p and compute 
S(l,L). Given the simple matrix product form of the ground-states, Eq. ([3]), it is 
possible to obtain explicit expressions for p using transfer matrix techniques. Though 
the reduced density matrix is large, it can be reduced to either a2x2or4x4 matrix, 
reflecting the fact that the allowed states for the subsystem A often is severely limited. 



3. Results 



3.1. S(l = 1,L), L even and odd 

We begin our calculations by looking at the case where /, the size of subsystem A, is 
one. We have found that of the four degenerate ground states, two (the S — 1, iSf. = ±1 
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Figure 1. The entanglement entropy of the system when I = 1, as a function of L, for 
the S = 1, Sfi = 1 state. The inset shows the convergence of the entropy towards its 
final value of C = —(2/3) ln(2/3) — (1/3) ln(l/3). The entropy converges exponentially 
with a constant of ln(3), which implies a correlation length of l/ln(3) in this system. 

states) have the same entanglement entropy by spin inversion symmetry. We therefore 
have three different cases to consider. 

3.1.1. S = 1, = ±1 In the following we explicitly consider the Sf. = 1 state of the 
triplet. We use Eqn. ([3]) to find the wave function and subsequently the reduced density 
matrix for each of the four ground states. We find that the reduced density matrix has 
only two eigenvalues, one of which is the probability, x, of a configuration containing 
the first spin (comprising all of subchain A) in a |+) state, and the remainder of the 
chain (subchain B) in a state with S^ B = 0, and the other is the probability, 1 — x, of 
the first spin being in a |0) state and the remainder of the chain having Sj. B = 1. For 
the = 1 state configurations with the first spin in a |— ) state does not occur at the 
AKLT point, hence, the probability, x is simply the on-site magnetization of the first 
element of the chain, (S[), which previously has been determined for both even and odd 
L:[33] 
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The reduced density matrix for the subsystem now takes the form: 





P 



x 




x 



This leads to the final equation for the Von Neumann entropy: 

S(l = 1,L) = -xln(x) - (1 -x)ln(l - x), (Any L) 



(7) 



(8) 



in agreement with the prior results by Alipour et al.[23] Interestingly, we have here 
related S(l = 1,L) directly to (Sf) an experimentally measurable quantity. The re- 
sult is a solution that converges exponentially fast with L towards a final value of 
C = —(2/3) ln(2/3) — (1/3) ln(l/3) as shown in Fig. [TJ The exponential form allows for 
a determination of a length scale which from Eqs. ([6]) and (jHJ) is seen to be £ = l/ln(3) 
equal to the known bulk correlation length of l/ln(3) [3Tj at the AKLT point. 



3.1.2. S = 0, = We now turn to a discussion of the first of the two S* = states. 
We note that the lower diagonal entry of Eq. (JHJ) is the spin inverse of the upper diagonal 
entry for even values of L. For odd L, the lower diagonal entry is the spin inverse of the 
upper diagonal entry times a factor of —1. This means that the singlet state when L is 
even, and the S% = triplet state when L is odd, are given by the trace of Eq. ([3]). In 
this case we find for I = 1 simply a constant independent of L: 

S(l = 1,L) = ln(3) (L even). (9) 

In this case the reduced density matrix is a 4x4 matrix with the following form: 



(10) 



Where in the case of I = 1, x is equal to 1/3. 

In the case where L is odd, the wave function is given by the upper left element 
of Eq. ([3]) minus the lower right element. This yields the following expression for the 
entanglement entropy of the singlet state when L is odd: 

(1 - (-3)- L+1 ^ 



X 














X 














X 














1 — X 



X 



(L odd) 



S(l 



3(1 -(-3)-, 
1, L) = -2x ln(x) - (1 - 2x) ln(l - 2x). 



In this case S(l = 1,L) now approaches the constant ln(3) exponentially with L. 



3.1.3. S — 1, iSy = Using the same arguments as for the 5 = case, we see that for 
even L the wave function is given by the upper left element of Eq. ()3]) minus the lower 
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right element, which produces the following result: 

3(W-V +L ) (LeV6n) 
S(l = 1, L) = -2x ln(ar) - (1 - 2x) ln(l - 2x), (12) 

the same result we found for odd L for the S = state. The reduced density matrix is 
still a 4x4 matrix, but it has the form: 



X 














X 














y 





_ 








1 - 2x 


= 0. 









y 



(13) 



Where for / = 1, y 

In the present case, for odd L, the wave function is given by the trace of Eq. ((3]), 
so the entanglement entropy is given by the formula of Hirano et al.: 

S(l — 1,L)— ln(3) (L odd) (14) 

To summarize, we have found that S(l = 1,L) approaches either —(2/3) ln(2/3) — 
(1/3) ln(l/3) or ln(3) in some cases in an exponential manner with L in other cases the 
result S(l = 1, L) is independent of L. 

3.2. S(l,L), L even 

We now generalize our results to any size I of the subsystem A, with L even. 

3.2.1. S — 1, iSj. — ±1 Since the entanglement is the same for Sj> = ±1 we in the 
following take Sj< = 1. In the case of Sj, = 1, we find that p again has two eigenvalues. 
One is the probability of finding the system in a state such that the total magnetization 
of subsystem A A = 1, and the total magnetization of subsystem B = 0. The 
other is the opposite case: S^ A = 0,S^ B = 1. If we denote the first eigenvalue by x, 
then the Von Neumann entropy is given by Eq. (jSJ), and the reduced density matrix has 
form of Eq. (J7J) The value x at any I is given by: 

;i - ( -3)-Q(i + ( -3)-^ 

2(1 -(-3)" 

S(l, L) = -x ln(x) - (1 - x) ln(l - x). (15) 

This result is plotted as a function of I when L = 12 in Fig. At L — 21, this gives 
a result of exactly S(l = L/2,L) = ln(2) independent of L. We also note that when 
I and L are both large, S(l, L) converges to ln(2) again in an exponential manner on 
a length scale of l/ln(3). This asymptotic value of the entanglement entropy seems 
natural since in the present case the boundary of the subsystem A will cut a single 
valence bond resulting in a contribution of ln(2) to the entanglement entropy. 



x = —1 , owm ( L eVen ) 
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Figure 2. The entanglement entropy for the S = 1, 5f. = 1 state as a function 



of / when L = 12 (Eq. (fT5|) ). The open circles indicate the corresponding exact 
diagonalization results used as a check. The inset shows the convergence of our data 
towards its final value of ln(2). The entropy therefore converges exponentially as 
/ -» L/2. 



3.2.2. S = 0, S T = In the case where L is even, the wave function is given by the 
trace of Eq. ([3]), so the reduced density matrix is given by Eq.(jT0l) and the entanglement 
entropy is given by the following formula: [H] 

(l-(-3)^)(l-(-3)" L+ 



x 



(L even) 



S(l,L) 



4(1 - (-3)-^) 
-3xln(x) — (1 — 3a;) ln(l — 3x) 



(16) 



Two of the degenerate eigenvalues correspond to the probabilities of subsystem A having 
a = an d the sum of the remaining degenerate eigenvalue and the non-degenerate 



eigenvalue is the probability of subsystem A having S TA = 0. Strikingly, in this case 
the entanglement entropy quickly approaches 2 ln(2) for large /, L, twice the result for 
the S = 1, S T = ±1 states. As above we can argue that this asymptotic value of the 
entanglement entropy is natural since in addition to cutting a single valence bond at the 
boundary, the subsystem now also cuts the singlet formed by the two effective S — 1/2 
chain boundary excitations resulting, in a contribution of 21n(2) to the entanglement 
entropy. It is possible to argue that the result, Eq. ffTB"]) . is independent of the boundary 
conditions and the above result does agree with previous results for periodic boundary 
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Figure 3. The entanglement entropy for the S = 0, Sj, = state (red) and the 
S = 1, Sj. = state (purple) as a function of I when L = 12 (Eq. (TT5)) ). (color online) 
The open symbols indicate the corresponding exact diagonalization results used as a 
check. The inset shows the convergence of the analytical data to the final value of 
21n(2). 



conditions. 

3.2.3. S = 1, Si? = In this case the reduced density matrix now has four eigenvalues 
and is given by Eq.(ITBl). The two degenerate eigenvalues correspond to the probability of 
subsystem A being in a state with Sj< A — ±1, The sum of the non-degenerate eigenvalues 
is the probability that S^ A — 0. We then find the following equation: 

(l-(-3)-Q(l-(-3)- L+t ) 

x = ; : — t- L even 

4(1 - -3 1 ^ 



-3)-')(l - (-3)- L+l+1 " 



V 



4(1 - (-3)"^) 
S(l,L) = -2xln(x) -yHy) 

-(l-2z-y)ln(l-2x-y). (17) 

Again we observe that the asymptotic value of the entanglement entropy is 21n(2). 
Hence, we have explicitly showed that the three states of the triplet S = 1, = 0, ±1 
only the two = ±1 states related by spin inversion yield the same entanglement 
entropy while the S = 1, = state not only differs by an overall factor of ln(2) but 



Exact Results for the Bipartite Entanglement Entropy of the AKLT spin-1 chain 



9 



also in subleading terms. The results Eqs. (TTtjj) and (TiTj) are plotted in Fig. [3] where 
they are compared with exact diagonalization results. Recent work [31] has suggested 
that the Haldane phase is characterized by two-fold degeneracy in the eigenvalues of the 
reduced density matrix of the ground state, which is consistent with these findings. 



3.3. S{l,L), L odd 

As before when we detailed the S{1 — 1, L) 
the parity of L for the general S(l, L). For 
the entanglement entropy also in this case. 



case we expect rather strong dependence on 
completeness we now give the equations for 



3.3.1. S — 1, jSj. = ±1 In this case there is no dependence on the parity of L and the 
result is the same as for even L given in Eq. ( fl5l) . 



3.3.2. S = 0, Sj. = By the same arguments as the / = 1 case, the entanglement 
entropy for odd L is given by: 

X ~ 4(1 -(-3)-*) ( d) 

(l-(-3)-Q(l-(-3)- £+i+1 ) 
y ~ 4(1 - (-3)"*) 

S(l,L) = -2xln(x) -y\n(y) 

-(l-2x-y)\n(l-2x-y). (18) 



3.3.3. S — 1, (Sy = Similarly, the wave function for S — 1, S?, — 0, L odd is given by 
the trace of Eq. (E]), and so the entanglement entropy is: 

4(1 - (-3)-^) (L ° dd) 

S(l, L) = -3x ln(x) - (1 - 3x) ln(l - 3z). (19) 



4. Conclusions 



We have obtained explicit analytical equations for the bipartite entanglement entropy 
of a spin-1 chain at the AKLT point for all four states of the ground-state manifold. For 
the case where S — 1, <Sy = ±1, we have found that for large system sizes the entangle- 
ment entropy approaches ln(2) while for the S = 0, 5f. = and S = 1, = cases, the 
entanglement entropy approaches 21n(2). Hence, the entanglement entropy is in this 
case not SU(2) invariant. In all cases where there is an explicit I or L dependence have 
we found that the asymptotic value is approached in an exponential manner defining a 
length scale of £ = 1/ ln(3) equal to the bulk correlation length. Also of interest are the 
time-reversal and spin-reversal invariant states \<f>) = \ + 1) ± | — 1). Exact calculations 
of the entanglement entropy of these states for small L suggests that their entanglement 
entropy converges towards 21n(2), which is consistent with the interpretation that each 
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factor of ln(2) corresponds to a cutting of a bond. We have so far been unable to obtain 
an explicit formula for the entanglement entropy of these states. 
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